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Abstract 

The (3+l)-dimensional Dirac equation with position dependent mass 
in 4-vector electromagnetic fields is considered. Using two over-simplified 
examples (the Dirac-Coulomb and Dirac-oscillator fields) , we report energy- 
levels crossing as a spectral property or as an effect of the hidden super- 
symmetric quantum mechanical language and/or quasi-parity signatures. 
Under different settings of the related interactions' way-of-coupling into 
Dirac equation, it is observed that the two ultimate/effective descendents, 
Dirac-Coulomb and Dirac-oscillator, exhibit different conditions on the 
energy- levels crossings. 

PACS numbers: 03.65.Ge, 03.65.Pm, 03.65. Fd, 03.65. Ca 

1 Introduction 

The search for exact-solvability for quantum mechanical systems, both non- 
relativistic (Schrodinger equation) and relativistic (e.g., Klein-Gordon and Dirac 
equations), is inviting and desirable. Exactly-solvable quantum mechanical sys- 
tems are vital ingredients for (among others that are mathematically motivated, 
say) the conceptual understanding of physics and for inspiring the structure 
of the numerical methods designed to solve more complicated physical prob- 
lems. However, whilst intensive attention was paid to exact-solvability of the 
non-relativistic Schrodinger equation, the relativistic Klein-Gordon and Dirac 
equations remained unfortunate and only partially attended in the yet already 
partially-explored Dirac territories. 

Some exactly-solvable potentials, for example, are known to belong to some 
distinctive classes of shape invariant potentials [1] . Within each class of which 
the so-called point-canonical-transformation (PCT) [2] would map the solution 
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(eigenvalues and cigcnfunctions) of one into another. On the other hand, su- 
persymmetric quantum mechanics [3] and potential algebras [4] (among others 
of course) are known to be used to obtain exact solutions for quantum mechan- 
ical systems. Yet, in between exactly-solvable and non-exactly-solvable there 
exists the gray zone of conditionally-exactly-solvable (i.e., all the spectrum is 
obtained) [5] and quasi-exactly-solvable (i.e., part of the spectrum is obtained) 
[6] potentials models. 

Nevertheless, we may recollect that the supersymmetric quantum mechanical 
language is realized as a hidden/built-in symmetry in the (l+l)-dimensional 
Dirac equation (cf., e.g. [3,7-9]). Nogami and Toyama [7] have reported that the 
associated Schrodingcr supersymmetric-partner Hamiltonians share the same 
energy spectrum including the lowest states unless Dirac equation allows a zero- 
mode (i.e., zero-energy bound-state). Jackiw and Rebbi [9] have, moreover, 
reported that only for some certain topological trends where the Lorentz scalar 
potential S(x) is localized (i.e., S(x) — > for x — > ±oo), Dirac equation would 
allow a zero-mode. 

One should be reminded, hereupon, that a Lorentz scalar potential (and/or 
an almost mathematically and partially-physically equivalent position-dependent 
mass) in Dirac equation is mainly motivated by the MIT bag model of quarks 
(cf., e.g., [10] and references therein). Yet, a Lorentz scalar potential couples 
to the mass of the fermion instead of its charge and the related positive and 
negative energies exhibit identical behaviors. Moreover, the position-dependent 
mass settings are useful models to study, for example, the energy density many- 
body problem, electronic properties of semiconductors and quantum dots, etc 
(cf., e.g., sample of references in [11-18]). We may also recollect that the phe- 
nomenon of energy-levels crossing is responsible for electron transfer in protein, 
it underlies stability analysis in mechanical engineering, and mathematically 
appears in algebraic geometry (cf., e.g., [19] and related references therein). 

Very recently [20], we have reemphasized the hidden/built-in supersymmet- 
ric quantum mechanical language in the spectrum of the (l+l)-dimensional 
Dirac equation with position-dependent mass and complexified Lorentz scalar 
interactions. We have reported the " quasi-parity" signature on the Dirac spec- 
trum and discussed energy- levels crossings related to supersymmetry and/or 
"quasi-parity" signatures. We have observed that the supersymmetric signa- 
ture on the (1+1)-Dirac spectrum is documented through the emergence of 
"exact" isospectral (i.e., including the lowest-state) partner Hamiltonians for 
" even" -quasi-parity, whereas the partner Hamiltonians share the same energy 
spectrum with a "missing" lowest-state for "odd" -quasi-parity, at least for the 
examples discussed therein. It would be interesting, we contemplate, if such 
studies are extended to cover the (3+l)-dimensional radial Dirac equation with 
different models of interactions' couplings and with position-dependent mass 
settings. Such studies merely exist in the literature, to the best of our knowl- 
edge, and may very well add a new flavour to the readily "multi-flavoured" Dirac 
equation. 

This article is organized as follows. In section 2, we recollect the (3+1)- 
dimcnsional radial Dirac equation with position-dependent mass in a four- vector 
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electromagnetic field. Therein, we realize that the decoupled one-dimensional 
Schrodinger-like radial Dirac equations exhibit supersymmetric language only 
when M (r) = V (r) = 0, where M (r) = m(r) + S (r) with the position- 
dependent mass m (r) , the Lorentz scalar potential S (r) , and the Lorentz vec- 
tor potential V (r). In section 3, the consequences of an equally- mixed non-zero 
Lorentz vector and Lorentz scalar potential settings are discussed through illus- 
trative examples: Dirac-Coulomb-I, Dirac-oscillator-I, and Dirac-oscillator-II. 
The consequences of V (r) ^ M (r) with the magnetic interaction A (r) = — 
C2 ( r ) / PC2 ( r )] 7^ ( sce Eq.(7) below) are given in section 4, along with il- 
lustrative examples: Dirac/Klcin-Gordon-Coulomb-II and Dirac/Klcin-Gordon- 
oscillator-III. In section 5, we report the consequences of an equally-mixed 
Lorentz vector and scalar "/ree "-fields (i.e., M (r) = V (r) = 0). The super- 
symmetric quantum mechanical language and the " quasi-parity" signatures on 
the spectra of a Dirac-oscillator-toy and a Dirac- Coulomb-toy models are re- 
ported in the same section. We conclude in section 6. 

2 Radial Dirac equation with position depen- 
dent mass in a 4-vector electromagnetic field, 
recollected 

The Hamiltonian describing a Dirac particle (in h = c = e = 1 units) in a 
four- vector electromagnetic field = (^A a , Aj = (Vi-Aj reads 

H D = a ■ (V - % ~T) + Pm + V, (1) 

with 

aj \ n ( 1 
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where aj are Pauli's 2x2 matrices and 1 is the 2x2 unit matrix. Under 
spherically symmetric settings, A — > f A (r) , V — > V (r) accompanied (for the 
convenience of the current study) by a position-dependent mass Lorentz scalar 
field, m — > m + m(r) + S (r) = m + M (r), the two-component Dirac equation 
reads 

/ m + M(r) + V(r) * + A (r) - 8 r \ ( g (r) \ _ ( g (r) \ 
{ f+A(r)+d r -m- M(r) + V(r) ) \ f (r) ) ~ \ f (r) ) 

(2) 

with energy E, and n in the centrifugal term is given by 

f -(^ + 1) forj = £+l/2 
k= { ' =► k(k+1)=£(£+1), (3) 

\ I for j =1-1/2 

where I = 0, 1, 2, • • • is the angular momentum quantum number. Equation (2) 
decouples into 

CAr)g(r)-(^+A(r)-d r ) / (r) = (4) 
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C 2 (r) f(r)-(^+A(r) + d r ) g(r) = 



where 



d ( r ) = (E - m) - V (r) - M (r) , 
C 2 (r) = (£ + m) - V (r) + M (r) . 
Substituting / (r) of (5) into (4) would, with 



A (r) = - + A (r) 
r 



imply 



+ i(r) 2 -I' (r) 



C 2 (r) 



i(r)+a r -C1WC2W 



)g(r) = 



(5) 

(6) 
(7) 

(8) 
(9) 



where primes denote derivatives with respect to r. Moreover, a substitution of 
the form 



g (r) = 2 (r) exp 



; i^(r) 



V (r) - AT (r) 
C 2 (r) 



(10) 



would remove the first order derivative and result in a one-dimensional Schrodingcr- 
like equation 



{-d 2 r + A {rf - A' (r) + U 2 (r) - d (r) C 2 (r)} </> 2 (r) = 



where 
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2C 2 (r)_ 



(11) 

(12) 



Similarly, substituting g (r) of (4) into (5) and taking 



/ ( r ) = 4>i (r) exp 
would imply 



Pi(r) 



; P[(r) = 



V (r) + M' (r) 



{-9 r 2 + i (r) 2 + i' (r) + C/r (r) - d (r) C 2 (r)} ^ (r) 



where 



4 Ui (r)>/ 2(iW 



(13) 
(14) 



It should be noted that the decoupled radial Dirac one-dimensional Schrodinger- 
likc equations, (11) and (13), possess natural (though hidden/built-in) super- 
symmetric quantum mechanical language only when M (r) = V (r) = 0. In 
the forthcoming experiments we shall be focusing on the upper component one- 
dimensional Schrodingcr-like radial Dirac equation in (11) and (12). 
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3 Consequences of an equally- mixed Lorentz vec- 
tor and scalar fields; V (r) = M (r) 

An equally-mixed Lorentz vector and scalar fields, V (r) = M (r) , would imply 
( 2 (r) = E + m, d (r) = E — m — 2M (r). Consequently, (11) reduces to 

{ -8 2 + A(r) 2 -A' (r) + 2M (r) [E + m] } 4> 2 (r) - [E 2 - m 2 ] ^ 2 (r) , (15) 
where A(r) = f + A (r), as given in (8). 

3.1 Dirac-Coulomb-I: 

For A (r) = a/r and M (r) = b/r equation (15) reads 

+ ^ + ^ (r) = [E 2 - m 2 ] ^ (r) , (16) 

where 



K = K ■ 



a + {j + 1/2) for k = + (j + 1/2) 
a - (j + 1/2) for K = - (j + 1/2) 



Equation (16) admits exact solution that can be very well inferred from the well 
known radial Schrodinger-Coulomb problem to yield (with the radial quantum 
number n r = 0, 1, 2, • • • ) 

[E 2 -m 2 ] =- b2[E + m] ; h = n r + K + l>0, (17) 



which in turn implies 

(n 2 - b 2 ) 
+ b 2 



lib \ IV V I , 

E = ( 18 ) 



However, the fact that k = o± (j + 1/2) would manifest energy- levels crossings 
to obtain. That is, a state labeled by hi = n r i + a + ji + 3/2 would cross with 
a state labeled by h 2 — n r2 + a — fa + 1/2 when 

h\ — b 2 h\ — b 2 _ 

hJTv = hJTv^ n2 = ni ^ nr2 - nrl=J1+j2 + 1 - 

Moreover, one may wish to mind the consequences associated with a com- 
plexified coupling constant in M (r) , and hence in V (r) , in such a way that 

V (r) — M (r) = b/r = —ib /r (e.g., simulating, say, the interaction of a point 
nucleus with an imaginary charge iZe and a particle of charge — e). In this case, 

V (r) = M (r) is "PT-symmetrized [21] and 

to (h 2 + b 2 ) 



E'PT 



bl 



Not only the energy spectrum E-pT follows similar energy-levels crossing sce- 
nario as that of (18), but also it suffers from the so called flown-away (cf., e.g., 
[21]) states that disappear from the spectrum when h = \b Q \ . 
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3.2 Dirac-Oscillator-I: 

For A (r) = br/2 and M (r) = equation (15) implies 

f *2 & 2 2 . M . ^ r^2 21 



</> 2 (r) , (19) 



where k (k + 1) = £ (£ + 1) for both k = - (£ + 1) ; j = £ + 1/2 and k = tj = 
I — 1/2 is considered. This would result in 

E 2 - m 2 = b(2n r + £ + 3/2) + nb- |, (20) 

to imply 

E± = ±yjm 2 + b (2n r + i + ~k + 1). (21) 

Obviously this result depends on the combination of the quantum numbers 
2n r + I = A and splits into 

± Vm 2 + b (A + j + 3/2) for k = + (j + 1/2) 
' = ^ ±Vm 2 +6(A-. 7 -l/2) for k = - (j + 1/2) (22) 

One should pay attention to the possible energy-levels crossings that occur 
between positive energy sets or negative energy sets. These are unavoidable 
energy-levels crossings manifested by k = ± (j + 1/2) and admits the following 
scenario: A state labeled by Ai and j\ crosses with a state labeled by A 2 and 
j2 for all b values when 

3 1 

Ai + ji + 2 = A 2 - h - 2 A 2 - A i = 3i + h + 2 

Nevertheless, quantum numbers related degeneracies are also feasible at different 
values of b. That is, when 

Ai = A 2 = A, ]\ ^ j 2 => A (b 2 - h) = &iji + 62.72 + 3bl ^ b2 
and when 

Ai ^ A 2 , ji = j 2 = j => A 2 &2 - A161 = j (h + b 2 ) + 3bl + b2 
3.3 Dirac-Oscillator-II: 

For A (r) = a/r and M (r) = B 2 r 2 /2 equation (15) reads 

{-9 2 r + + B 2 [E + m] r 2 } ^ 2 (r) = [E 2 to 2 ] ^ 2 (r) , (23) 

In this case 

[E 2 - m 2 ] = 2BVE + to (2n r + k + 3/2) (24) 
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which would lead to E = — m (to be discarded) and, with N± = 2n r + a ± 
(j + 1/2) + 3/2 > 0, 

£ = -™+Qd /3 + 2m£ ± 1/3 ) 2 (25) 

where 



e± = 27B iV± + 3^/ -24m3 + 81B 2 J\£; N ± > ^ ^ (26) 

In this case, it is obvious that energy-levels crossings occur when £ + = One 
may, for the sake of simplicity, choose the case where 



, 8m 3 

which would imply that a state labeled by n r \ and j\ crosses with a state labeled 
by n r2 and j 2 , for all a, when 

£+ = t_=^N+= N_=^ n r2 - n rl = (ji + j 2 + 1) /2 

4 Consequences of A (r) = — C 2 (r) / [2£ 2 (r)] 7^ 
and V (r) 7^ M (r) 

If we consider the class of interactions where A (r) = — C 2 ( r ) /2C2 ( r ) 7^ an d 
V(r) 7^ M (r), Dirac equation in (11) reduces to 

-d 2 r + ^H±A - Cl (r) C 2 (r) U 2 (r) = 0, (27) 



which is, in fact, in exact form as that of the radial one-dimensional Klein- 
Gordon (KG) equation 

-31 + l l^L -[E-V (r)] 2 + [m + M (r)A <f> 2 (r) = 0. (28) 



Moreover, it should be noted hereby that the case where A (r) = —n/r— 
(2 (r) /2£ 2 (r) 7^ and V (r) 7 M (r) corresponds to the s- waves (i.e., I = 0) 
solution of (28). Hence, one need not consider it as a separate case to deal with. 

4.1 Dirac/Klein-Gordon-Coulomb-II: 

For V(r) = a\/r and M (r) = a 2 /r equation (28) reads 

+ + 2Kg r +Q2W] } 2 (r) = [E> - ™ 2 ] 2 (r) , (29) 
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with 

£ = - l - + ^j(l+l/2) 2 -a 2 +a 2 > 0, 

and admits exact solution of the form 

r-^o 9i \ot\E + a 2 m] 2 , , 
[E 2 - m 2 ] = - 1 ' J ; N = n r + C + 1 > 0. 

This would lead to 

E± 



m 



M 2 + a \ 



± 



U\Ol 2 



N 2 -a 2 



-1 1/2 



N 2 + a 2 ) Af 2 + a 2 



to yield, for various especial cases of coupling constants, 

± [l-al/N 2 ] 1 ' 2 \ N>\a 2 \ for Qi=OA« 2 ^0 

= < ± [l + a{/N 2 Y 112 for ai^0A« 2 = 

(-a 2 ± TV 2 ) / (TV 2 + a 2 ) for a x = a 2 = a ^ 



E± 



m 



(30) 



(31) 



(32) 



It should be obvious that energy-levels crossings for (32) are not feasible at all. 
At this point, one should notice that the effect of positive/negative k is absent 
in the process of choosing A (r) = — ( 2 (r) /2( 2 (r) ^ and V (r) ^ M (r). 
However, the complexification of the coupling constants would manifest flown- 
away states (cf., e.g., [21] ) to obtain and disappear from the spectrum. For 
example, for the case 

a± = a 2 = ia ^ y = 



in 



M 2 



the energy states fly-away and disappear from the spectrum when J\f = \a \, 
and for 




12 ' 



energy states fly-away when J\f = \a[\. However, flown-away states never occur 
for the case 

E± IN 2 - a 12 

ai - 0, a 2 = ia 2 =► — = ±y ^ 2 , 

but rather the system loses its observability (i.e., reality and discreteness) and 
collapses when \a' 2 \ > N. Nevertheless, in connection with such complexified 
version of a\ and a 2 , the reader may consult Mustafa [21] for more comprehen- 
sive details on the spectral properties of a similar complexified Coulombic fields 
(but with 7 in [21] replacing C in the current example). 
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4.2 Dirac/Klein-Gordon-Oscillator-III: 

For V (r) = 0, and M (r) = (3 1 /r+ [3 2 r — m equation (28) implies 



-8 2 + V 2 ' + f3 2 2 r 2 + 2/3^2 \ cj> 2 (r) = E 2 2 (r) 




where 




(33) 



with 



(34) 



Obviously, neither energy-levels crossings nor flown-away states are feasible for 
this model. Nevertheless, one should pay attention to the parametric settings 
that may lead to imaginary energies and consequently system collapse. 



5 Consequences of equally-mixed "/ree"-fields: 
supersymmetry and quasi-parity 

It is obvious that with equally- mixed Lorentz vector and Lorentz scalar "free"- 
fields (i.e., V (r) = M (r) = 0) one would obtain, from (11) and (13), 



stands for an effective supersymmetric-like partner potentials. The supersym- 
metric language is very well pronounced in this process, therefore. Yet, the 
quasi-parity shall emerge in the forthcoming two Dirac-toy models. 

5.1 A Dirac-oscillator-toy 

Let us consider a "toy" model: 



that results in an effective supersymmetric partner "Dirac-oscillator-toy" po- 
tential of the form 



where A± = E 2 - 



{-d 2 + V ± (r)}0 ± {r) = X±<f> ± (r), 
-m 2 , (f> + (r) = cj) 1 (r), (f>_ (r) = <j> 2 (r), and 
V± (r) = V u (r) = A (r f + luA' (r) ; uj = ±1 



(35) 



A 1 

A (r) = + -Br; R3i,B>0, 

r 2 




(36) 



9 



In the repulsive/ attractive-like core, moreover, one may replace A (A + u>) by 
a (a + 1). In this case, a would denote quasi-angular momentum quantum num- 
ber and a = —1,0 may very well correspond to "even" and "odd" quasi-parity 
(i.e., q = (— 1) <J+1 ), respectively. For more details on quasi-parity convention 
the reader may refer to, e.g., Mustafa and Znojil [22] and Znojil [23] and related 
references cited therein. In a straightforward manner, however, one can show 
that 

a(a + l)=A(A + u)=>a = ~+q{A+^j. (37) 

Under such settings, it is obvious that both supcrsymmetric-like partner poten- 
tials in (36) admit exact closed form solutions (cf., e.g., Mustafa and Znojil [22], 
and Znojil [23]): 

A w = -(4n r + 2qA + u)q + 2)- b(a~ |) ; n r = 0, 1, 2, • • • , (38) 
which would split (with u = ±1) into 

2B (n r + 1) 
2B (n r -A+ if^ 

2Bn r 

2B (n r -A+ ±)(40) 

We may now pay attention to the supersymmetric language "signature" which 
is documented in the facts that \+. q= -i = A_ i?= _i, for odd quasi-parity, and 
A + g=+ i = A_ :(J=+ i + const. = A_ i?= +i + 2_B, for even quasi-parity. That is, 
for an even quasi-parity the superpartner potentials possess identical spectra 
with a missing lowest state, whereas for an odd quasi-parity the superpartner 
potentials are "exactly" isospectral. Nevertheless, the consequence of such "/wd- 
den"-supersymmetry in Dirac equation is very well pronounced in the related 
Dirac spectra (with A± = E 2 — to 2 ): 

E +tQ=+1 = + ^/m 2 + 2B (n r + l) 

< / — o / '. TV ; n r = 0, 1,2, • • • (41) 

£+,,=-! = +^m 2 + 2B(n r - A+\) ' V ' 

{E-. q=+1 = -^m 2 + 2Bn r 
E_^ = -^m 2 + 2B(n r -A + l) ^ = °' ^ ' ' ' ^ 

In this case, the quasi-parity signature appears in the energy-level crossings 
between the two sets of energies in (41) and between those in (42). That is, the 
two sets of energies in (41) cross with each other when 

E + (n r = n rl ,q = +1) = E + (n r = n r2 ,q = -1) => n r2 - n rl = A + i (43) 



22 

\+, q = — (4n r + 2qA + q + 3) - BA 



A_, 9 = —(An r + 2qA-q + l)-BA 




E +,q - 

E-, g = 
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and those in (42) cross with each other when 



(n r = n r s, q = +1) = £_ (n r = n r4 , q = -1) n r4 - n r3 = A - -. (44) 

5.2 A Dirac-Coulomb-toy 

On the other hand, a "toy" model of the form 

Air) = + B; R 9 A,B > 0, 

r 

would result in an effective supersymmetric partner "Dirac-Coulomb-toy" po- 
tentials 



which admits exact solution of the form 

\2 



Au — — 



(ABY 



77/ 



+ B ; ft = n r + a + 1 > 0. n r = 0, 1, 2, • • 



(45) 



(46) 



This result would split into 



A+, 9 = B 2 <1 — A 2 



B <1 — A 



n r + q I A 



n r + q[A-- 



1 -2 



(47) 
(48) 



Each of which, respectively, yields 



^+,9 



9=+l 



{ 



S 2 1 - A 2 



A- 



I]" 2 } 



X+, q =-i = B 2 {l-A 2 [n r -A]- 2 } 

' X-, q =+i = B 2 {l- A 2 [n r + A}- 2 } 
\-, q =-i = B 2 {l - A 2 [n r - A + If 2 } 



(49) 



(50) 



Yet, it should be noted here that both even and odd quasi-parity eigenvalues 
A+ i9= _i and A_ ig= +i allow zero-modes (i.e., zero-energy) at n r = and there- 
fore they do not share the same spectrum (a Nogami's and Toyama's [7] ob- 
servation). However, A+ !(J= +i and A_ !(J=+ i do not allow zero-modes and hence 
they have identical spectra but with a missing lowest state in A_ >g=+ i. Similar 
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trend is also observed for A+. 9= _i and A_ ig= _i where the lowest state is missed 
in A+ j9= _i. Consequently, the corresponding Dirac spectra read 



E. 



E+,q=+l 



= + Jm 2 + £? 2 jl - A 2 [n r + A + If 2 } 



£+, 9 =-i = + Jm 2 + B 2 {l - A 2 [n r - A}' 2 } 



(51) 



= < 



= -\ m 



B 2 {l-^ 2 K + Af 2 } 



[l-A 2 ln r -A+l]- 2 } 



(52) 



Evidently, energy- levels crossings obtain between the two sets of energies in (51) 
or between the two sets of energies in (52). That is, a state E +iq=+ \ (n r = n r i) 
crosses with a state £7 +i9= _i (n r = n r2 ) when 

n r i + A + 1 = n r2 — A => n r2 — n r \ = 2A + 1, 

and E_ )g=+ i (n r = n r3 ) crosses with (n r = n r4 ) when 

n r 3 + A = n r 4 — A + 1 n r 4 — n r 3 = 2 A — 1. 

Nevertheless, the energy sets in (51) and (52) loose their reality and become 
pure imaginary in the following manner: 

E +tq=+1 e Cform 2 |B 2 < A 2 B 2 / (n r + A + I) 2 , 

E+. q= ^ e C for m 2 + B 2 < A 2 B 2 / (n r - A) 2 , 

£_,„=+! £ C for m 2 + B 2 < A 2 B 2 / (n r + A) 2 , 

S_, 9= _! e Cform 2 + B 2 < A 2 B 2 / {n r - A + I) 2 . 



6 Conclusion 

The inspiration of the current work is stimulated by our subsequent study [20] 
of the hidden/built-in supersymmetric quantum mechanical language and/or 
quasi-parity signatures on the spectrum of the (1+1)-Dirac equation. However, 
as long as Dirac and Klein-Gordon wave equations are concerned, the energy- 
levels crossing phenomenon/paradox (the discussion of which already lies far 
beyond our current proposal, cf., e.g., [19] for more details) as a spectral property 
or as a consequence of the supersymmetric language and/or quasi-parity is left 
an almost-forgotten one. Our purpose, even with the current over-simplified 
radial Dirac/Klein-Gordon examples, was to fill this gap at least partially. 

In the light of the current study, we have observed that under different 
settings of the magnetic interaction field A (r) , and likewise the related inter- 
actions' way-of-coupling into Dirac equation, the two ultimate/effective descen- 
dents, Dirac-Coulomb and Dirac-oscillator, perform energy-levels crossing at 
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different conditions. Moreover, our observations in section 5 on the spectral 
properties of the radial supersymmetric partner Hamiltonians (i.e., A± i9= ±i in 
(39), (40), (49), and (50)) re-confirm Nogami's and Toyama's [7] ones on the 
(l+l)-dimcnsional Dirac equation. 

Finally, we contemplate that the variety of settings (presented in this work) 
of the related interactions' way-of-coupling (i.e., Eqs. (15), (27), and (35)) 
into the (3+l)-dimensional radial Dirac equation would enrich the number of 
exactly/quasi-exactly/conditionally-exactly solvable Dirac models. Yet, the so- 
lution of the most general radial case in (13) still resides in the mathemati- 
cally challenging Hcrmitian, non-Hcrmitian, and pscudo-Hcrmitian [20,21,24-26] 
Dirac territories. 
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